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ANALYSIS  OF  AEROSOL  AGING  IN  THE  ROTATING  DRUM  CHAMBER 


1.  INTRODUCTION 

The  rotating  drum  aerosol  chamber  was  developed  to  maintain  aerosol 
concentrations  over  long  periods  of  time  in  order  to  study  aerosol  aging.  By  rotating  the  drum 
about  a  horizontal  cylindrical  axis  we  can  avoid  gravitational  settling  which  removes  significant 
fractions  of  the  aerosol  from  static  aerosol  chambers  over  times  that  are  large  compared  to  the 
chamber  height  divided  by  the  Stokes  settling  velocity.  The  rotating  drum  is  useful  for  studying 
the  viability  of  biological  aerosols  over  long  times  and  for  studying  forces  and  effects  that 
normally  would  be  overwhelmed  by  gravitational  settling.  Two  theoretical  treatments  of  aerosol 
mechanics  in  a  rotating  drum,  (Gruel  et  al.  1987)  and  (Asgharian  and  Moss  1992),  will  be 
compared  to  one  another  and  referred  to  as  Gruel’s  theory  and  Moss’s  theory.  Then  predictions 
from  these  theories  will  be  compared  to  recent  rotating  drum  chamber  measurements  (Sutton 
2003).  Attempts  will  be  made  to  explain  any  significant  differences  between  measurements  and 
predictions  using  these  theories. 


2.  EXPERIMENTAL  DETAILS 

In  the  experimental  study  (Sutton  2003)  aluminum  oxide  aerosol  was  introduced 
into  a  1  meter  diameter  by  1  meter  long  cylindrical  chamber  (drum)  from  a  large  plastic  trash  bag 
containing  aerosol  disseminated  pneumatically  using  the  SRI  sonic  nozzle.  This  was  done, 
starting  with  nearly  zero  initial  chamber  volume,  by  moving  the  drum  piston  wall  away  from  the 
opposite  wall  creating  negative  pressure  that  pulled  the  aerosol  from  the  bag  into  the  drum 
through  a  section  of  Tygon  tubing  attached  to  the  lA  inch  diameter  port  located  at  the  center  of 
the  opposite  drum  wall  on  the  drum  axis.  Aerosol  samples  were  taken  to  determine  aerosol 
concentration  and  particle  size  as  a  function  of  time,  rotation  rate  and  sampling  port  location  by 
moving  the  drum  piston  wall  toward  the  opposite  wall  creating  positive  pressure.  Small  volumes 
of  aerosol  were  pushed  through  one  of  four  Vi  inch  diameter  sampling  ports,  in  the  opposite  wall 
located  at  radial  positions  between  the  axis  and  5  cm  from  the  curved  wall,  into  an  Aerodynamic 
Particle  Sizer  manufactured  by  TSI.  Essentially  the  chamber  was  filled  and  emptied  like  a 
syringe. 


3.  THEORETICAL  BASIS 

Gruel’s  theory  describes  aerosol  particle  spiral  trajectories  in  the  laboratory  frame 
of  reference  due  to  combined  gravitational  and  centrifugal  forces.  A  suspension  volume  is 
defined  for  each  particle  size  as  the  fraction  of  aerosol  remaining  after  one  revolution  of  the 
chamber  has  cleared  away  all  aerosol  within  a  distance  of  one  radius  of  gyration  of  the  particle. 
This  is  best  understood  in  the  rotating  drum  frame  of  reference  where  particles  gyrate  due  to 
rotating  gravitational  force  and  intercept  the  wall  if  within  a  distance  of  one  radius  of  gyration. 
As  they  gyrate  due  to  gravity  they  slowly  move  radially  out  toward  the  drum  wall  due  to 
centrifugal  force  following  spiral  paths  in  the  laboratory  frame  of  reference  and  flat  helical  paths 
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with  small  pitch  angles  due  to  the  relatively  small  centrifugal  force  compared  to  gravitational 
force  in  the  rotating  drum  frame  of  reference.  Gruel’s  theory  is  appropriate  when  centrifugal 
force  is  much  less  than  gravitational  force  predicting  aerosol  total  number  suspended  after  one 
revolution  and  decay  thereafter  as  aerosol  particles  move  into  the  clear  zone  due  to  centrifugal 
force. 


Moss’s  theory  uses  the  same  equations  of  motion  but  introduces  the  concept  of 
limiting  trajectory  to  describe  the  combined  effects  of  gravity  and  centrifugation  over  the 
complete  range  of  rotation  rates.  Moss’s  theory  predicts  aerosol  concentration  decay  beginning 
when  drum  rotation  starts. 


4.  RESULTS  AND  ANALYSIS 

The  radius  of  gyration  of  a  particle  in  the  rotating  drum  chamber  is 


(1.1) 


where  r,  g  and  co  are  respectively  the  particle  relaxation  time,  acceleration  due  to  gravity  and 
angular  velocity  of  drum  rotation.  The  relaxation  time  can  be  expressed 


-  -  PdlCc 
18  TJ 


(1.2) 


in  terms  of  77,  p,  d  and  Cc  which  are  respectively  the  dynamic  viscosity  of  air,  the  particle 
density,  diameter  and  Cunningham  slip  correction  factor 


C=\  +  Kn 


1.257 +  0.4  exp 


LH 

V  Kn 


(1.3) 


where  Kn  =  2/1 /d  is  the  Knudsen  number  which  is  twice  the  ratio  of  the  mean  free  path  of  air 

molecules,  0.071  pm  at  standard  temperature  and  pressure,  over  the  particle  diameter.  The  total 
number  of  aerosol  particles  suspended  in  the  chamber  decays  exponentially  according  to  Gruel’s 
theory  so  that  the  fraction  of  aerosol  remaining  is 


N_ 

N0 


exp  (-2^7) 


(1.4) 


where  R  and  t  represent  the  drum  radius  and  time  beginning  after  one  revolution.  Equivalent 
expressions  for  the  fraction  of  aerosol  remaining  based  on  Moss’s  theory  describe  what  happens 


10 


during  the  first  revolution  and  include  the  possibility  that  centrifugal  acceleration  exceeds 
gravitational  acceleration 


-^-  =  exp(-2 Tco2t^j  for  -^-<l-exp (-r&r/) 


TV  l  +  exp(-2r6Tf)  j 


Nn 


n 


iL  i_[^l 


.  H, 

+  sin  — L  + 
R 


H  ^  h  (H-X 

0-[T)  + 


f  H  \ 

exp(-2rni2t)sin'1 


for 


' 

(  n  \ 

H  ] 

1-exp 

-TOO  t 

1  +  exp 

-T(02t 

(1.5) 


—  =  0  for  —  >  l  +  exp(-r<y2t) 


where 


H  =  —J  1-2  cos(©r)  exp(-r<w2/)  +  exp(-2m>2?)  for  cot  <n 

co 

H  =  — [l  +  exp(-zw2/)]  for  mt>n 


Hx 


R2  +H2  -  R2  exp(-2rco2t) 


2  H 


H2=H-Hl 


(1.6) 


and  t  here  represents  the  time  beginning  when  the  first  revolution  begins. 

The  fraction  of  aerosol  remaining  suspended  in  the  drum  based  on  Moss’s  theory 
is  surface  plotted  in  Figures  1  through  13.  The  radius  of  gyration  and  fraction  of  aerosol 
remaining  after  one  revolution  based  on  Gruel’s  theory  is  surface  plotted  in  Figures  14  and  15. 
The  ratio  of  Gruel’s  solution  divided  by  that  of  Moss  predicting  total  aerosol  mass  or  number  is 
surface  plotted  in  Figures  16  through  28  after  compensating  for  the  1/rpm  time  offset  by 
assuming  aerosol  is  removed  at  a  constant  rate  during  the  first  revolution  to  create  a  clear  zone. 

It  should  be  pointed  out  that  these  expressions  for  total  number  of  aerosol  particles  suspended  in 
the  chamber  can  be  used  to  predict  total  aerosol  mass  suspended  and  to  predict  aerosol 
concentration  that  remains  homogeneous.  If  the  theoretical  predictions  of  total  number  of 
aerosol  particles  in  the  chamber  are  to  be  compared  to  concentration  measurements  taken  from  a 
single  port  location,  the  concentration  profile  within  the  chamber  must  be  known. 

The  rotating  drum  used  in  the  experiment  is  shown  in  Figure  29  and  was  loaned  to 
us  by  Owen  Moss.  The  drum  has  a  750  liter  capacity  and  four  sampling  ports  located  in  the 
drum  wall  opposite  the  piston  wall  are  shown  in  Figure  30.  Port  1  is  located  5  cm  from  the  outer 
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curved  wall  of  the  chamber  and  ports  2,  3,  and  4  are  located  14,  24,  and  50  cm  from  that  wall 
with  port  4  on  the  drum  axis.  Aerosol  concentration  spatial  distribution  can  be  monitored  by 
comparing  data  taken  from  these  four  ports  and  Figure  31  shows  that  the  initial  concentration 
and  size  distribution  are  homogeneous.  Figures  32  through  36  summarize  these  measurements 
for  1  pm  through  5  pm  aerodynamic  diameter  particles  respectively.  In  each  of  these  figures  the 
data  are  grouped  according  to  rotation  rates  of  0,  1,  2,  4,  6,  and  8  rpm  indicated  along  the  x  axis 
at  the  time  of  measurement.  The  ratios  of  number  concentrations  measured  through  ports  1,  2, 
and  3  over  those  measured  through  port  4  (the  axis  port)  are  shown  on  the  y  axis  and  define  the 
error  bars  at  times  t  =  0,  2,  4,  and  6  hours.  Concentrations  appear  to  be  fairly  homogeneous 
except  for  the  6  rpm  test  after  4  hours  with  axis  port  concentration  slightly  higher  than  the  other 
three  ports.  This  is  to  be  expected  considering  that  positive  pressure  in  the  drum  forces  the 
aerosol  through  the  port  from  radial  positions  surrounding  the  port  so  that  the  existence  of  a  clear 
zone,  even  at  port  1  near  the  wall,  will  not  be  evidenced  by  anything  like  a  step  function  change 
in  concentration. 

Concentration  fraction  measurements  taken  through  the  drum  axis  port  are  plotted 
at  all  rates  of  rotation  as  “percent  of  initial  particles”  in  Figures  37  through  42  for  respectively 
1, 2,  3, 4,  5,  and  6  pm  aerodynamic  diameter  aerosol  particles.  Corresponding  plots  based  on 
Moss’s  theory  predictions  are  shown  in  Figures  43  through  48.  Significant  differences  between 
predictions  and  measurements  are  obvious  especially  at  rpm=0  where  Moss’s  theory  greatly 
underestimates  suspension  times.  The  stirred  settling  model  is  appropriate  not  only  for  turbulent 
mixing  due  to  fans  but  also  convective  mixing  due  to  small  differences  between  air  and  chamber 
wall  temperatures  generally  encountered  with  static  chambers  (Fuchs  1964).  The  stirred  settling 
model  assumes  that  convective  or  turbulent  diffusion  maintains  uniform  concentration 
throughout  the  chamber  volume  up  to  the  laminar  boundary  layer  near  the  wall.  Within  the 
laminar  boundary  layer  diffusion  is  damped  and  Stokes  settling  becomes  the  dominate  mode  of 
particle  transport  to  the  wall  where  it  is  deposited  reducing  total  aerosol  particle  number  in 
accordance  with  the  equation 


N 

^=exn- 


2v/^i 

nR 


(1.7) 


for  a  tube  or  cylindrical  chamber  lying  on  its  side  where  v  =  rg  is  the  Stokes  settling  velocity. 
Measurements  at  rpm  =  0  are  summarized  in  Figure  49  and  the  stirred  settling  model  predictions 
appearing  in  Figure  50  are  a  good  match. 

Measured  aerosol  decay  rates  are  dramatically  reduced  when  the  drum  is  rotated 
suggesting  that  convective  diffusion  is  reduced  because  transport  near  the  wall  due  to 
gravitational  gyration  is  not  much  different  from  that  due  to  gravitational  settling  in  a  static 
chamber.  Convective  motion  of  buoyant  air  in  the  chamber  can  be  expected  to  involve  gyrations 
that  drive  progressively  smaller  eddies  until  viscosity  damps  out  the  process  and  drum  rotation 
will  smooth  out  temperature  gradients  in  the  drum  wall.  If  convective  diffusion  is  the  dominant 
mechanism  controlling  transport,  then  we  can  estimate  what  happens  by  applying  the  solution  for 
diffusion  inside  a  spherical  chamber  toward  a  perfectly  absorbing  walls  and  obtain 
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(1.8) 


N  6  ^exp (Dn2j2t/R2) 
\T  ~  ^2  2-1 


j= 1 


J 


where  D  and  R  are  respectively  the  convective  diffusion  coefficient  and  chamber  radius  (Fuchs 
1963).  Convective  diffusion  coefficients  can  be  expected  to  decrease  approaching  the  wall  until 
they  are  damped  out  at  the  laminar  boundary  layer  making  this  solution  irrelevant  for  a  static 
chamber.  However  when  the  drum  is  rotated  and  convective  diffusion  is  reduced  sufficiently  by 
buoyant  air  gyration  and  rapid  viscous  dissipation  of  small  eddies,  then  gravitational  gyration  of 
aerosol  particles  can  clear  a  zone  extending  one  radius  of  gyration  away  from  the  wall.  Because 
this  clear  zone  extends  well  beyond  the  laminar  boundary  layer  we  can  anticipate  the  possibility 
that  convective  diffusion  might  dominate  transport  by  centrifugal  drift  into  the  clear  zone  where 
gravitational  gyration  becomes  the  dominant  transport  mechanism.  This  can  be  incorporated  into 
the  diffusion  solution  simply  by  reducing  chamber  radius  by  the  radius  of  gyration 


N  _  6  ^e*p{Dx2j2t/(R-r0)2) 

~  2-J 


N, 


71 


j= 1 


J 


(1.9) 


calculations  using  this  equation  appear  in  Figure  51  for  the  largest  and  smallest  radii  of  gyration 
corresponding  respectively  to  d  =  6  pm  at  1  rpm  and  d  =  1  pm  at  8  rpm  for  convective  diffusion 
coefficients  of  1 0-2 , 1 0  3  and  1 0  4  cm2 /s  giving  concentrations  roughly  in  the  range  measured. 
These  six  curves  show  little  dependence  on  particle  size  and  rotation  rate  resulting  in  essentially 
three  curves  defined  by  the  three  diffusion  coefficients  which  are  intrinsically  independent  of 
particle  size  over  the  range  of  interest  (Fuchs  1964)  and  probably  independent  of  rotation  rate. 
Therefore,  the  observed  particle  size  and  rpm  dependence  in  concentration  decay  rate  is  not 
explained  by  the  diffusion  model.  These  convective  diffusion  coefficients  are  smaller  than  those 
encountered  in  the  stirred  settling  model  representing  the  rpm  =  0  case  where  diffusion  maintains 
essentially  homogeneous  concentration  in  spite  of  gravitational  settling. 

When  transport  into  the  clear  zone  by  diffusion  and  drift  due  to  centrifugal  force 
are  comparable  in  magnitude,  we  look  to  the  steady  state  one  dimensional  solution  (Friedlander 
1977)  for  insight.  The  deposition  flux  at  the  clear  zone  boundary,  where  n  =  nb,  is 


J=- 


VM. 


l-exp[-vc(/?-r0)/  D] 


(1.10) 


the  ratio  J(D  =  0)  /  J  of  deposition  flux  due  to  centrifugal  drift  alone  over  that  due  to  combined 
diffusion  and  drift  is  plotted  in  Figures  52  and  53  for  convective  diffusion  coefficients  of 
0.1  and  0.01  cm2 /s  over  a  range  of  particle  sizes  and  rates  of  rotation.  It  is  evident  that  diffusion 
which  is  intrinsically  independent  of  particle  size  will  have  the  greatest  influence  on 
progressively  smaller  particle  transport  as  the  rate  of  rotation  increases  such  that 
vc[R-r0)/D  <  1.  While  convective  diffusion  significantly  speeds  up  the  concentration  decay 
rate  of  smaller  particles  at  lower  rates  of  rotation  and  will  help  predict  the  measurements  for 
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these  particles,  it  has  a  negligible  effect  on  the  larger  particles  especially  at  the  higher  rates  of 
rotation  and  will  not  help  reconcile  the  large  differences  between  Moss’s  theory  predictions  and 
the  measurements. 

Another  pertinent  result  for  combined  drift  plus  diffusion  is  found  in  the  transient 
solution  for  aerosol  transport  through  an  electrostatic  precipitator  (Shapiro  1998)  in  the  limit 
where  the  Peclet  number,  drift  velocity  multiplied  by  transport  distance  divided  by  turbulent 
diffusion  coefficient  is  greater  than  10.  Here  diffusion  slows  the  deposition  process  according  to 
the  equation 


n 

—  =  exp 


vV 

ad) 


(1.11) 


and  works  to  decrease  rather  than  increase  the  aerosol  concentration  decay  rate  further  increasing 
differences  between  measurements  and  Moss’s  theory  predictions. 

Although  aerosol  particle  number  concentrations  were  on  the  order  of  -700/  cm’ 
during  the  drum  experiments,  the  long  residence  times  in  the  chamber  could  allow  coagulation  to 
play  a  role  along  with  centrifugal  drift  in  reducing  aerosol  concentration.  To  develop  this 
approach  we  note  that  homogeneous  concentration  is  maintained  outside  the  clear  zone  in  the 
rotating  drum  regardless  of  convective  diffusion  if  initial  concentration  is  homogeneous.  This  is 
because  the  rate  of  change  of  aerosol  concentration,  equal  to  the  rate  of  change  of  aerosol  mass 
in  a  volume  element  (flux  in  minus  flux  out)  divided  by  that  volume  element,  is  independent  of 
radius  for  centrifugal  acceleration  which  is  proportional  to  radius.  The  total  number  of  aerosol 
particles  in  the  rotating  drum  can  be  written 


N  =  n(R-r<j)2  Ln  (1.12) 

where  L  and  n  are  the  drum  length  and  homogeneous  aerosol  number  concentration.  The  rate  of 
change  in  total  number  is  equal  to  the  flux  of  particles  entering  the  clear  zone 

«-2  z(R-r„)Lvcn  (1.13) 

dt 

where  vc  =  co2  (i?  -  r0 )  r  is  the  particle  velocity  due  to  centrifugal  force.  Combining  the  two 
previous  equations  we  get 


dn  2v 

—  =  7 — ~^:n 
dt  ( R-r0 ) 


(1.14) 


or  simply 
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(1.15) 


dn 

~dt 


2co2Tn 


The  combined  effects  of  loss  due  to  centrifugal  drift  which  is  proportional  to  the 
number  concentration  and  coagulation  which  is  proportional  to  the  number  concentration 
squared  are  represented  by  the  differential  equation 

-^-  =  pn  +  Kn1  (1.16) 


for  a  homogeneous  cloud  which  has  the  solution  (Fuchs  1 964) 


1 


1  K 
—  +  — 


\no 


exp(fit)- 


£ 

P 


(1.17) 


where  n0  and  K  are  the  initial  concentration  and  coagulation  constant  with  /3  =  la/r .  In  the 

limit  where  coagulation  becomes  negligible,  K  — »  0  this  expression  provides  the  solution  for 
aerosol  concentration  as  a  function  of  time  when  loss  is  due  to  centrifugal  drift  into  the  clear 
zone.  Over  the  particle  size  range  of  interest  Brownian  coagulation  sets  a  lower  limit  of 
~10  9  cm  js  and  differential  gravitational  settling  (in  our  case  gyration)  sets  an  upper  limit  of 
~10  8  an'  js  for  the  coagulation  constant  or  more  precisely  for  the  collision  frequency  function 
(Friedlander  1977).  The  ratio  of  the  number  concentration  remaining  after  48  hours  with  over 
without  coagulation  is  surface  plotted  in  Figures  54  and  55  as  a  function  of  particle  size  and  rpm 
for  the  lower  and  upper  limit  values  of  the  coagulation  constant.  The  influence  of  coagulation 
can  be  seen  to  be  significant  in  reducing  number  concentration  of  smaller  particles  at  all  rates  of 
rotation  and  larger  particles  at  the  slower  rates  of  rotation  but  has  negligible  effect  on  the  number 
concentration  of  larger  particles  at  higher  rotation  rates.  Therefore,  differences  between 
measurements  and  Moss’s  theory  predictions  for  larger  particles  at  higher  rotation  rates  are  not 
explained. 


During  these  experiments  the  SRI  nozzle  disseminated  aluminum  oxide  powder 
by  aspirating  the  powder  into  and  through  an  aluminum  tube  using  a  sonic  Velocity  air  flow 
surrounding  the  tube  outlet  and  moving  along  the  tube  axis  away  from  the  tube.  This  same 
airflow  then  projects  the  powder  out  in  a  jet  that  dilutes  the  aerosol  concentration  while  inducing 
turbulence  that  deagglomerates  aerosol  particles.  The  aluminum  tube  is  grounded  and  it  is  likely 
that  triboelectric  charging  of  the  powder  occurred  as  the  powder  moved  through  the  tube  and 
scraped  the  tube  wall.  A  monopolar  charged  cloud  produced  in  this  way  would  expand  because 
of  mutual  repulsion  and  aerosol  concentration  would  decay  due  to  a  combination  of  centrifugal 
and  electrostatic  drift  into  the  zone  cleared  by  gravitational  gyration.  The  component  of  aerosol 
particle  velocity  at  the  edge  of  the  clear  zone  due  to  electrostatic  repulsion  between  that  particle 
and  the  monopolar  charged  cloud  is 


15 


(1.18) 


where  Fjm  is  the  acceleration  due  to  the  electrostatic  force,  Fe ,  acting  on  an  aerosol  particle  of 
mass  m  .  The  electrostatic  force  is 


(1.19) 


where  q  is  the  charge  on  the  particle  and  Q  is  the  total  charge  of  a  spherical  monopolar  charged 
cloud.  If  we  assume  that  the  amount  of  triboelectric  charge  is  proportional  to  the  surface  area  of 
the  particle,  then 


q  =  %ed2  (1.20) 

where  %d2  is  the  number  of  elementary  charges  e  =  4.8  T10'10es«  on  the  particle  in  the  cgs  unit 
system.  Assuming  homogeneous  concentration,  the  total  charge  contained  in  the  cloud  is 

Q  =  ~x(R-rof  n£e(d2)  (1.21) 

where  (d1^  reflects  an  average  surface  area  per  particle  which  depends  on  the  instantaneous  size 

distribution.  The  diameter  corresponding  to  yj(d2^  is  called  the  diameter  of  the  average  surface 
and  can  be  related  to  the  more  common  number  median  diameter 

dnmj  =  )  exp  (-In2  ag)  ,  (1 .22) 

surface  median  diameter 

=^)  exp  (in2  crg)  (1.23) 

and  mass  median  diameter 

dmmd=$f)™  P(21n2crg)  (1.24) 

descriptors  (Reist  1993)  for  a  log  normal  size  distribution  where  a g  is  the  geometric  standard 
deviation. 
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The  equation  controlling  aerosol  concentration  when  coagulation,  drift  due  to 
centripetal  acceleration  and  drift  due  to  triboelectric  charging  occur  simultaneously  in  a 
homogeneous  cloud  becomes 


dn  ..  2  2v 

- -  Bn  +  Kn  +- - h 

dt  H  (A-r0) 


(1.25) 


which  can  be  written 


-—  =  j3n  +  K'n2  (1.26) 

dt 


with  the  solution 


n  = 


1 


\no 


K ...  K 

H - exp  (Bt) - 

P )  P 


(1.27) 


where 


K  =K+— 
n 


(1.28) 


and 


2ve  S7TT^2d2  (d2^ e2n 
{R-r0)  3  m 


(1.29) 


Figure  56  shows  the  ratio  of  aerosol  remaining  after  24  hours  with  triboelectric 
charging  of  10 d2  electron  charges  per  particle  (d  in  microns)  over  that  remaining  without 
triboelectric  charging.  Initial  concentration  is  700  particles  per  cubic  centimeter,  (d2^  =  9 

microns  and  coagulation  is  ignored.  The  magnitude  of  for  this  case  is  plotted  in  Figure  57 

and  exceeds  K  over  almost  the  entire  size  distribution  range.  Figures  58  through  63  show  the 
combined  effects  of  centrifugal  and  electrostatic  drift  at  this  charging  level  into  the  clear  zone 
along  with  coagulation  at  the  Brownian  level  K  =  10~9  cm3 /s  assuming  the  coulombic  repulsive 
force  reduces  coagulation  to  that  level.  If  the  van  der  Waals  forces  holding  the  particles  together 
after  they  make  contact  do  not  exceed  the  repulsive  force  then  the  coagulation  constant  would 
have  to  be  set  to  zero  but  it  has  little  effect  on  the  results  at  the  Brownian  level.  Because  (d2  ^ 

continually  decreases  as  the  charged  cloud  decays,  we  found  that  by  simply  setting  it  equal  to  d2 
the  measurements  could  be  matched  better  than  by  using  any  fixed  value.  Thus  the  smaller 
particles  which  are  suspended  longer  on  average  encounter  a  repulsive  cloud  of  particles  minus 
the  larger  particles.  The  larger  particles  suspended  for  shorter  times  see  a  cloud  of  all  particle 
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sizes  but  the  larger  particles  have  much  greater  charge  per  particle  and  represent  the  majority  of 
charge  existing  in  the  cloud.  Predictions  using  this  approach  are  in  reasonable  agreement  with 
the  measurements  and  a  closer  fit  would  no  doubt  be  possible  by  further  adjusting  the  charging 

coefficient  £  or  (d2^ . 


5.  CONCLUSIONS 

The  stirred  settling  model  was  found  suitable  for  the  drum  chamber  when  it  is  not 
rotated.  Small  temperature  differences  between  the  wall  and  aerosol  drive  convective  diffusion 
that  is  sufficient  to  uniformly  mix  particle  up  to  10  pm  aerodynamic  diameter  even  in  a  one 
meter  high  chamber.  Once  the  drum  is  rotated  at  rates  of  1  to  8  rpm  convective  diffusion  is 
damped  out  by  convective  flow  gyration  and  viscous  dissipation  of  small  eddies.  Transport  is 
then  governed  by  centrifugal  drift  into  a  clear  zone  as  described  by  Gruel’s  and  Moss’s  theories. 
This  zone  extends  one  radius  of  gyration  away  from  the  wall  is  rapidly  cleared  by  gravitational 
gyration  in  one  revolution.  During  the  experiment  electrostatic  forces  due  to  triboelectric 
charging  within  our  dissemination  device  were  found  to  be  important  in  determining  the  aerosol 
decay  rate  at  the  rates  of  rotation  studied.  Triboelectric  charging  in  the  dissemination  device 
used  has  always  led  to  negligible  effects  during  static  aerosol  chamber  tests.  This  is  because  the 
electrical  forces  involved  are  small  compared  to  gravity  in  controlling  transport  through  the 
laminar  boundary  layer  and  small  compared  to  convective  or  turbulent  diffusion  controlling 
transport  outside  the  laminar  boundary  layer.  The  rotating  drum  essentially  replaces  gravity  with 
centrifugal  force  and  at  the  same  time  damps  out  convective  diffusion.  Instead  of  gravitational 
settling  through  the  laminar  boundary  layer  the  rotating  drum  has  gravitational  gyration  through 
the  laminar  boundary  layer  and  a  clear  zone  extending  one  radius  of  gyration  away  from  the  wall 
that  is  created  in  one  revolution.  In  order  to  study  other  phenomena  involving  forces  that  are 
small  compared  to  gravity,  triboelectric  charging  should  be  avoided  by  neutralizing  the  aerosol 
as  it  exits  the  dissemination  device  or  the  dissemination  method  should  be  altered  to  avoid 
triboelectric  charging. 
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Figure  2.  Aerosol  Rotating  Drum  (R  =  50  cm,  t  =  24  hr) 
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Figure  4.  Aerosol  Rotating  Drum  (R  =  50cm,  rpm  =  .5) 
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50 

Figure  9.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  10g) 


Figure  10.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  5g) 


23 


Figure  11.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  2g) 
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Figure  13.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  ,5g) 
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Figure  20.  Aerosol  Rotating  Drum  (R 


Figure  21 .  Aerosol  Rotating  Drum  (R  =  50  cn 


Figure  22.  Aerosol  Rotating  Drum  (R 
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Figure  23.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  20g) 


Figure  24.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  10g) 
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Figure  25.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  5g) 
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Figure  27.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  1  g) 


50 


Figure  28.  Aerosol  Rotating  Drum  (R  =  50  cm,  D  =  .5g) 
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Figure  32.  Ratio  of  Measured  Concentration  of  1  gm  Particles  in  Ports  1,  2,  and  3 

Over  that  in  Port  4  (Axis  Port)  at  Various  Times  and  as  a  Function  of  RPM 
Rotation  Rate  Plotted  on  the  Abcissa. 


Figure  33.  Ratio  of  Measured  Concentration  of  2  gm  Particles  in  Ports  1,  2,  and  3 

Over  that  in  Port  4  (Axis  Port)  at  Various  Times  and  as  a  Function  of  RPM 
Rotation  Rate  Plotted  on  the  Abcissa. 
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Figure  34.  Ratio  of  Measured  Concentration  of  3  gm  Particles  in  Ports  1 ,  2,  and  3 

Over  that  in  Port  4  (Axis  Port)  at  Various  Times  and  as  a  Function  of  RPM 
Rotation  Rate  Plotted  on  the  Abcissa. 


Figure  35.  Ratio  of  Measured  Concentration  of  4  gm  Particles  in  Ports  1,  2,  and  3 

Over  that  in  Port  4  (Axis  Port)  at  Various  Times  and  as  a  Function  of  RPM 
Rotation  Rate  Plotted  on  the  Abcissa. 
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Figure  36.  Ratio  of  Measured  Concentration  of  5  |im  Particles  in  Ports  1,  2,  and  3 

Over  that  in  Port  4  (Axis  Port)  at  Various  Times  and  as  a  Function  of  RPM 
Rotation  Rate  Plotted  on  the  Abcissa. 


One  Micron  Diameter 


Figure  37.  Concentration  of  1  gm  Diameter  Particles  Measured  through  Axis  Port 


Two  Micron  Diameter 


Figure  38.  Concentration  of  2  gm  Diameter  Particles  Measured  through  Axis  Port 
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Figure  39.  Concentration  of  3  jam  Diameter  Particles  Measured  through  Axis  Port 


Figure  40.  Concentration  of  4  gm  Diameter  Particles  Measured  through  Axis  Port 


Figure  41.  Concentration  of  5  gm  Diameter  Particles  Measured  through  Axis  Port 
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Six  Micron 


Figure  42.  Concentration  of  6  gm  Diameter  Particles  Measured  through  Axis  Port 


N/No  Hoss  Theory  (rpm  =  0  ,1,2,4,6,8;D=1^) 


Figure  43.  Asgharian  and  Moss  Theory  Predictions  for  1  gm  Particles 


N/No  Hoag  Theory  ( rpm  =  0 ,  i ,  2 , 4  ,  6 , 8 ;  D  =  2#0 


Figure  44.  Asgharian  and  Moss  Theory  Predictions  for  2  gm  Particles 
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H/No  Moss  Theory  (  rpm.  =0 , 1 2 , 4,  6 ,8  ;  D=3^j 


Figure  45.  Asgharian  and  Moss  Theory  Predictions  for  3  gm  Particles 


N/Nc  Hubs  Theui*  (  rpm  =  0  ,  .  ,  2 , 4 , 6  ,  8  ;  D  =  4 ft ) 


Figure  46.  Asgharian  and  Moss  Theory  Predictions  for  4  gm  Particles 


N/No  Moss  Theory  ( rpm  =  0 , 1 , 2 , 4 , 6  ,  8  ;  D  =  5u) 


Figure  47.  Asgharian  and  Moss  Theory  Predictions  for  5  gm  Particles 
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N/No  Hoss  Theory  (  rpa  =  0 , 1 ,  2 , 4 , 6  ,  8  ;  D  =  6  ju) 


Figure  48.  Asgharian  and  Moss  Theory  Predictions  for  6  pm  Particles 


Summary  of  RPM=0  Data 


0  12  3  4  5  6  7 

Time  (hr) 


Figure  49.  Summary  of  Static  (rpm  =  0)  Drum  Measurements  through  Axis  Port 
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Figure  50.  Stirred  Settling  Theory  Predictions  of  1, 2,  3,  4,  and  5  pm  Diameter 

Aerosol  Concentration  Fractions  Remaining  as  a  Function  of  Time  in  Hours 


N  No 


Figure  51.  Diffusion  Theory  Predictions  of  Concentration  Fraction  Remaining  as  a  Function 
of  Time  in  Hours.  Three  Pairs  of  Curves  Corresponding  to  Diffusion  Coefficients 
of  0.01,  0.001,  and  0.0001  cm2/s  at  the  Largest  and  Smallest  Radii  of  Gyration 
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*5 


1 


N/N< 


Figure  61.  Ncharge  =  10*DA2  (rpm  =  1,  2,  4,  6,  8;  D  =  4p) 


N/No 


Figure  62.  Ncharge  =  10»DA2  (rpm  =  1,  2,  4,  6,  8;  D  =  5|li) 


N/No 


Figure  63.  Ncharge  =  10*DA2  (rpm  =  1,  2,  4,  6,  8;  D  =  6p) 
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APPENDIX 


MATHEMATICA  CODE  GENERATING  SURFACE  PLOTS 


MATHEMATICA  file  MossN(D,RPM).nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*. 07 1  *  1 0A-4/d)( 1 .257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  IfTPi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fi  =p; 

f2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

0  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
thour  =  3; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  f2]],  {DD,  0.01, 

20},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,t=3hr)", 

AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "N/No"},  PlotPoints  ->  20]; 
thour  =  24; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  f2]],  {DD,  0.01, 

20},  {rpm,  .0001,2}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50em,t=24hr)", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 
thour  =  96; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {DD,  0.01, 

20},  {rpm,  .0001,  2}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,t=96hr)", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 
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gl  =  Plot3D[h/r,  {DD,  0.01,  20},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,t=3hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "h/r"},  PlotPoints  ->  20]; 


MATHEMATICA  file  MossN(T,D).nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour* 3 600; 

CC  =  1  +  (2*  .07 1  *  1 0A-4/d)(  1.257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  lf}Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

f2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5 
p*ArcSin[h2/r]); 

f3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm=  .5; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If} hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001, 

48},  {DD,  .01,20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=5)", 
AxesLabel  ->  {"thour",  "D(p)",  "N/No"},  PlotPoints  ->  20]; 
gl  =  Plot3D[h/r,  {thour,  0.001,  48},  {DD,  .01,  20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=.5)", 
AxesLabel  ->  {"thour",  "D(p)M,  "h/r"},  PlotPoints  ->  20]; 
rpm  =  1; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001, 

48},  {DD,  .01,20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm— 1)", 
AxesLabel  ->  {"thour",  "D(p)",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  Iffhrhigh  <=  h/r,  D,  f2]],  {thour,  0.001, 
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48},  {DD,  .01,20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=2)", 
AxesLabel  ->  {"thour",  "D(p)",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  5; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  (thour,  0.001, 
48},  (DD,  .01,20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=5)", 
AxesLabel  ->  ("thour",  "D(p)",  "N/No"},  PlotPoints  ->  20]; 


MATHEMATICA  file  MossN(T,RPM).nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fi  =  p; 

f2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5 
p*ArcSin[h2/r]); 

f3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +pA0.5; 

DD  =  20; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  Iffhrhigh  <=  h/r,  f3,  f2]],  (thour,  0.0001, 

50},  (rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=20p)", 
AxesLabel  ->  ("thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  10; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  fi,  f2]],  (thour,  0.0001, 

50},  (rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=10p)", 
AxesLabel  ->  ("thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  5; 
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gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.0001, 

50},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=5p)", 

AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  2; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.0001, 

50},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=2p)", 

AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  1; 
gl  -  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.0001, 

50},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=lp)", 

AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  .5; 
gl  -  Plot3D[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.0001, 

50},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=.5p)", 
AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 
gl  =  Plot3D[h/r,  {thour,  0.0001,  50},  {rpm,  .0001,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=.5p)", 
AxesLabel  ->  {"thour",  "rpm",  "h/r"},  PlotPoints  ->  30]; 

MATHEMATICA  file  MossD=l.nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC=  1  +(2*.071*10A-4/d)(1.257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 

pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r=  50; 
fl  =p; 

f2  -  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5 
p*ArcSin[h2/r]); 
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£3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm  =  .0000001; 

DD  =  1; 

gO  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g2  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g3  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g4  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g5  =  Plot[If{h/r  <=  hrlow,  fl,  If{hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Show[g0,  gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {0, 1}, 

PlotLabel  ->  "  Rotating  Drum  (rpm=0,l,2,4,6,8;D=l\[Mu])"]; 

MATHEMATICA  file  MossD=2.nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour* 3 600; 

CC  =  1  +  (2  *  .07 1  *  1 0A-4/d)( 1.257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(1 8*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
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£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

O  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm  =  . 0000001; 

DD  =  2; 

gO  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=2\[Mu])", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm=  1; 

g2  =  Plot[If[h/r  <=  hrlow,  fl,  Ifflirhigh  <=  h/r,  G,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=2\[Mu])", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g3  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  G,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D:=2\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm=:4; 

g4  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  G,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=2\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If{hrhigh  <=  h/r,  G,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=2\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g5  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  G,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=2\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Show[g0,  gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Rotating  Drum  (rpm=0,l,2,4,6,8;D=2\[Mu])"]; 


MATHEMATICA  file  MossD=3.nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*. 07 1  *  1 0A-4/d)(  1.257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 
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omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +p)/2- 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

0  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm  =  .0000001; 

DD  =  3; 

gO  =  Plot[lf[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g2  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g3  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  fi,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g4  =  Plot[If{h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g5  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=3\[Mu])M, 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Show[g0,  gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Rotating  Drum  (rpm=0, 1 ,2,4,6,8;D=3\[Mu])"] ; 


MATHEMATICA  file  MossD=4.nb 

rho  =  1 ; 

mu=  1.83*1 0A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
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t  =  thour*3600; 

CC  =  1  +  (2*. 07 1  *  1 0A-4/d)(l  .257  +  .4*Exp[-l.  l*d/(2*.  071  *10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

0-0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm  =  .0000001; 

DD  =  4; 

gO  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  £2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=4\[Mu])", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm=  1; 

g2  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  f2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm-l,D=4\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"} ,  PlotPoints  ->  20]; 
rpm  =  2; 

g3  =  Plot[If[h/r  <=  hrlow,  fl,  If{hrhigh  <=  h/r,  B,  £2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=4\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g4  -  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  B,  f2]],  {thour,  0.001, 25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=4\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  B,  f2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=4\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g5  =  Plot[If[h/r  <-  hrlow,  fl,  If[hrhigh  <=  h/r,  B,  £2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=4\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Show[g0,  gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Rotating  Drum  (rpm=0, 1 ,2,4,6,8;D=4\[Mu])"]; 
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MATHEMATICA  file  MossD=5.nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*. 07 1  *  1 0A-4/d)(  1 .257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

f3  —  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm  =  .0000001; 

DD  =  5; 

gO  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  f2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D“5\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g2  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  (thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=5\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g3  -  Plot[lffh/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  O,  £2]],  (thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm:=2,D=5\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g4  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  (thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=5\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  (thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=5\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g5  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  £3,  £2]],  (thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=5\[Mu])", 
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AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Show[g0,  gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Rotating  Drum  (rpm=0,l,2,4,6,8;D=5\[Mu])"]; 


MATHEMATICA  file  MossD=6.nb 

rho  =  1 ; 

mu=  1.83*1 0A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*  .07 1  *  1 0A-4/d)(  1.257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

£3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rpm  =  .0000001; 

DD  -  6; 

gO  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001,  7}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g2  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001,  7}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g3  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  {thour,  0.001,  7}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g4  =  Plot[If[h/r  <=  hrlow,  fl,  Iffhrhigh  <=  h/r,  £3,  £2]],  {thour,  0.001,  7}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  {thour,  0.001, 7}, 
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PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g5  =  Plot[If[h/r  <=  hrlow,  fl,  Ifphrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001,  7}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

ShowfgO,  gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {0,  1), 

PlotLabel  ->  "  Rotating  Drum  (rpm=0,l,2,4,6,8;D=6\[Mu])"]; 


MATHEMATICA  file  RatioGruelMoss(T,RPM).nb 


rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour* 3 600; 

CC  =  1  +  (2*. 07 1  *  1 0A-4/d)(  1 .257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  *  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +p)/2- 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

f3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh=  1  +pA0.5; 
rO  =  tau  980/omega; 
fgl  =  If[thour  <=  1/60/rpm, 

1  -  (1  -  (r  -  r0)A2/rA2)(thour/60/rpm)Exp[-2*tau*t*omegaA2],  (r  - 
rO)A2  Exp[-2*tau*t*omegaA2]/rA2]; 

fg2  =  0; 

fg  =  If[r  <=  rO,  fg2,  fgl]; 

DD  =  20; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/f3,  fg/f2]],  {thour, 

0.0001,  50},  {rpm,  .1,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=20|i)", 

AxesLabel  ->  {"thour",  "rpm",  "Ng/Nm"},  PlotPoints  ->  30]; 

DD  =  10; 
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gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/B,  fg/G]],  {thour, 
0.0001,  50},  {rpm,  .1,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=10p)", 
AxesLabel  ->  {"thour",  "rpm",  "Ng/Nm"},  PlotPoints  ->  30]; 
DD  =  5; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/B,  fg/f2]],  {thour, 
0.0001,  50},  {rpm,  .1,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=5|j.)", 
AxesLabel  ->  {"thour",  "rpm",  "Ng/Nm"},  PlotPoints  ->  30]; 
DD  =  2; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/B,  fg/G]],  {thour, 
0.0001,  50},  {rpm,  .1,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=2|i)", 
AxesLabel  ->  {"thour",  "rpm",  "Ng/Nm"},  PlotPoints  ->  30]; 
DD  =  1; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/B,  fg/G]],  {thour, 
0.0001,  50},  {rpm,  .1,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=lp)", 
AxesLabel  ->  {"thour",  "rpm",  "Ng/Nm"},  PlotPoints  ->  30]; 
DD  =  .5; 
gl  =Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/B,  fg/G]],  {thour, 
0.0001,  50},  {rpm,  .1,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,D=.5p)", 
AxesLabel  ->  {"thour",  "rpm",  "Ng/Nm"},  PlotPoints  ->  30]; 


MATHEMATICA  file  RatioGruelMoss(T,D).nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*  .07 1  *  1 0A-4/d)( 1 .257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 

pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2H!h); 
h2  =  h-hl; 
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r  =  50; 
fl  =p; 

£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

f3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +pA0.5; 
rO  =  tau  980/omega;  fgl  = 

If[thour<=  1/60/rpm, 

1  -  (1  -  (r  -  r0)A2/rA2)(thour/60/rpm)Exp[-2*tau*t*omegaA2],  (r  - 
rO)  A2  Exp  [-2  *  tau*  t*  omegaA2]/rA2] ; 
fg2  =  0;  ' 

fg  =  If[r<=rO,  fg2,  fgl]; 
rpm  =  .5; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/f3,  fg/f2]],  {thour,  0.001, 
48},  {DD,  .01,  20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=.5)", 

AxesLabel  ->  {"thour",  "D(p)",  "Ng/Nm"},  PlotPoints  ->  20]; 
rpm  =  1; 
gl  =  Plot3D[ 

If{h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/f3,  fg/f2]],  {thour,  0.001, 
48},  {DD,. 01,  20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=l)", 

AxesLabel  ->  {"thour",  "D(p)",  "Ng/Nm"},  PlotPoints  ->  20]; 
rpm  =  2; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/fi,  fg/f2]],  {thour,  0.001, 
48},  {DD,. 01,  20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=2)", 

AxesLabel  ->  {"thour",  "D(p)",  "Ng/Nm"},  PlotPoints  ->  20]; 
rpm  =  5; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/f3,  fg/f2]],  {thour,  0.001, 
48},  {DD,  .01,20}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,rpm=5)", 

AxesLabel  ->  {"thour",  "D(|x)",  "Ng/Nm"},  PlotPoints  ->  20]; 


MATHEMATICA  file  RatioGruelMoss(D,RPM).nb 

rho  =  1 ; 

mu=  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
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t  =  thour*3600; 

CC  =  1  +  (2*.07 1 *  1 0A-4/d)(l  .257  +  .4*Exp[- 1 . 1  *d/(2*.07 1  *  1 0A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  IftPi  <= 

omega*t,  (tau*980/omega)(l  + 
pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A.5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

f3  =  0; 

hrlow  =  1  -  pA0.5; 
hrhigh  =  1  +  pA0.5; 
rO  =  tau  980/omega; 
fgl  =  If[thour  <=  1/60/rpm, 

1  -  (1  -  (r  -  r0)A2/rA2)(thour/60/rpm)Exp[-2*tau*t*omegaA2],  (r  - 
rO)  A2  Exp  [-2  *  tau*  t*  omega  A2]/r A2  ] ; 

fg2  =  0; 

fg  =  If[r  <=  rO,  fg2,  fgl]; 
thour  =  3; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/O,  fg/f2]],  {DD,  0.01, 
20},  {rpm,  .01,  5}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,t=3hr)", 

AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "Ng/Nm"},  PlotPoints  ->  20]; 
thour  =  24; 
gl  =  Plot3D[ 

If[h/r  <=  hrlow,  fg/fl,  Ifthrhigh  <=  h/r,  fg/fi,  fg/f2]],  {DD,  0.01, 
20},  {rpm,  .01,2}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,t=24hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "Ng/Nm"},  PlotPoints  ->  20]; 
thour  =  96; 
gl  =  Plot3D[ 

If{h/r  <=  hrlow,  fg/fl,  If[hrhigh  <=  h/r,  fg/f3,  fg/f2]],  {DD,  0.01, 
20},  {rpm,  .01,2}, 

PlotLabel  ->  "Aerosol  Rotating  Drum  (R=50cm,t=96hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "Ng/Nm"},  PlotPoints  ->  20]; 


rho  =  1; 


MATHEMATICA  file  GruelN(T,D)&RPMopt.nb 
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g  =  980; 

mu  =  1.83*10A-4; 
d  -  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  -  thour* 3 600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
rO  =  tau  980/omega; 
fi  =  Iff  thour  <=  1/60/rpm, 

1  -  (1  -  (r  -  r0)A2/rA2)(thour/60/rpm)Exp[-2*tau*t*omegaA2],  (r  - 
rO)A2  Expf-2*tau*t*omegaA2]/rA2]; 

£2  =  0; 

omegaopt  =  (g/2/r/t)A0. 333333  +  tau*g/3/r; 
rpmopt  =  60*omegaopt/2/Pi; 
r  =  50; 
rpm  =  0.5; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {DD,  .001,  20},  {thour,  0.0001,  50}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,rpm=0.5)", 
AxesLabel  ->  {"D(p)M,  "thour",  "N/No  "},  PlotPoints  ->  20]; 
rpm  =  1; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {DD,  .001,  20},  {thour,  0.0001,  50}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,rpm=l)", 
AxesLabel  ->  {"D(p)",  "thour",  "N/No  "},  PlotPoints  ->  20]; 
rpm  =  2; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {DD,  .001,  20},  {thour,  0.0001,  50}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,rpm=2)", 
AxesLabel  ->  {"D(p)",  "thour",  "N/No  "},  PlotPoints  ->  20]; 
gl  =  Plot3D[rpmopt,  {DD,  .001,  20},  {thour,  0.0001,  50}, 

PlotLabel  ->  "Optimum  Rotation  Rate  (R=50cm)", 

AxesLabel  ->  {"D(p)",  "thour",  "RPMopt  "},  PlotPoints  ->  20]; 


MATHEMATICA  file  GruelN(D,RPM).nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  ,4*Expf-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
rO  =  tau  980/omega; 
fl  =  Iff  thour  <=  1/60/rpm, 

1  -  (1  -  (r  -  r0)A2/rA2)(thour/60/rpm)Exp[-2*tau*t*omegaA2],  (r  - 
rO)A2  Exp[-2*tau*t*omegaA2]/rA2]; 


APPENDIX 


59 


£2  =  0; 
r=  50; 
thour =  3; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {DD,  0.001,  20},  (rpm,  .01,5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,t=3hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 
thour  =  24; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  (DD,  0.001, 20},  {rpm,  .01,5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,t=24hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 
thour  =  48; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {DD,  0.001, 20},  {rpm,  .01,5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,t=48hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 


MATHEMATICA  file  GruelN(T,RPM).nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 

tau  =  rho*dA2  CC/(18*mu); 

p  =  Exp[-2*tau*t*omegaA2]; 

rO  =  tau  980/omega; 

fl  =  If[  thour  <=  1/60/rpm, 

1  -  (1  -  (r  -  r0)A2/rA2)(thour/60/rpm)Exp[-2*tau*t*omegaA2],  (r  - 
rO)A2  Exp[-2*tau*t*omegaA2]/rA2]; 

£2  =  0; 

r  =  50; 

DD  =  1; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {thour,  .001,  50},  {rpm,  .01,  5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,D=lp)", 

AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  2; 

gl  =  Plot3D[If{r  <=  rO,  £2,  fl],  {thour,  .001,  50},  {rpm,  .01,  5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,D=2p)", 

AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  5; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {thour,  .001,  50},  {rpm,  .01,  5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,D=5p)", 

AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 

DD  =  10; 

gl  =  Plot3D[If[r  <=  rO,  £2,  fl],  {thour,  .001,  50},  {rpm,  .01,  5}, 
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PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,D=10p)", 
AxesLabel  ->  {"thour",  "rpm",  "N/No"},  PlotPoints  ->  30]; 
DD  =  20; 

gl  =  Plot3D[If[r  <=  rO,  f2,  fl],  (thour,  .001,  50},  (rpm,  .01,  5}, 
PlotLabel  ->  "Gruel  Rotating  Drum  (R=50cm,D=20p)", 
AxesLabel  ->  ("thour",  "rpm",  "N/No"},  PlotPoints  ->  30, 
PlotRange  ->  (0,  1}]; 


MATHEMATICA  file  StirredSettling.nb 

rho  =  1; 

mu  =  1.83*10A-4; 

d  =  DD*10A-4; 

t  =  thour* 3 600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 

tau  =  rho*dA2  CC/(18*mu); 

f  =  Exp[-2*tau*980*t/r/Pi]; 

r  =  50; 

DD  =  1; 

gO  =  Plot[f,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=l\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 

DD  =  2; 

gl  =  Plot[f,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 

DD  =  3; 

g2  =  Plot[f,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=3\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 

DD  =  4; 

g3  =  Plot[f,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 

DD  =  5; 

g4  =  Plot[f,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=5\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 

Show[g0,  gl,  g2,  g3,  g4,  PlotRange  ->  ((0,  25},  (0,  1}}, 

PlotLabel  ->  "  Drum  Stirred  Settling  (D=l,2,3,4,5\[Mu])"]; 

Show[g0,  gl,  g2,  g3,  g4,  PlotRange  ->  {(0,  7},  (0,  1}}, 

PlotLabel  ->  "  Drum  Stirred  Settling  (D=l,2,3,4,5\[Mu])"]; 


MATHEMATICA  file  ElectrostaticDispersion.nb 

th=  l/(4*pi*B*nuA2*eeA2*nO); 
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B  =  l/(3*pi*eta*d); 
ee  =  4.8*10A-10; 
eta  =  1.83*10A-4; 
d  =  dmicron*10A-4; 
nO=  1000; 
thhour  =  th/3600; 
nu  =  nufact*dmicronA2; 

gl  =  Plot3D [thhour,  {nufact,  1,  10},  {dmicron,  1,  6}, 

PlotLabel  ->  "Charged  Cloud  Halflife(n=  1 000", 

AxesLabel  ->  {"\[Nu]",  ”D(\[Mu])",  "t(hr)  "},  PlotPoints  ->  20]; 
nu  =  nufact*dmicron;  gl  = 

Plot3D [thhour,  {nufact,  5,  15},  {dmicron,  1,6}, 

PlotLabel  ->  "Charged  Cloud  Halflife(n=1000", 

AxesLabel  ->  {"\[Nu]",  "D(\[Mu])",  "t(hr) "},  PlotPoints  ->  20]; 
nu  =  nufact;  gl  = 

Plot3D[thhour,  {nufact,  10,  30},  {dmicron,  1,6}, 

PlotLabel  ->  "Charged  Cloud  Halflife(n=1000", 

AxesLabel  ->  {"\[Nu]",  "D(\[Mu])",  "t(hr)  "},  PlotPoints  ->  20]; 

MATHEMATICA  file  ElectrostaticDispersionFinal.nb 

rho  =  1; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  .4*Exp[-l.l*d/(2*.  071  *10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
h  =  If[Pi  <= 

omega*t,  (tau*980/omega)(l  + 

pA.5),  (tau*980/omega)(l  -  2*Cos[omega*t]*pA.5  +  p)A5]; 
hi  =  (r*r  +  h*h  -  r*r*p)/(2*h); 
h2  =  h  -  hi; 
r  =  50; 
fl  =p; 

£2  =  (1  +  p)/2  - 

l/Pi*(hl/r*(l  -  (hl/r)A2)A.5  +  ArcSin[hl/r]  +  h2/r*(p  -  (h2/r)A2)A.5  + 
p*ArcSin[h2/r]); 

£3  =  0; 

hrlow  =  1  -  pA0.5; 

hrhigh  =  1  +  pA0.5; 

nn  =  l/(Exp[2*tau*t*omegaA2]* 

(l/nn0  +  KK/(2*tau*omegaA2))  -  KK/(2*tau*omegaA2)); 

KK  =  6*10A-10  +  8*Pi*tau*qqA2/3/mm; 

Feoverm  =  KK*nn; 
mm  =  rho*Pi*dA3/6; 
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qq  =  ncharge*4.8*10A-10; 
ncharge  =  10*DDA2; 
nnO  = 1000; 
ff  =  nn/nnO; 
rpm  =  .0000001; 

DD  =  1; 

gl  =  Plot[If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g2  =  Plot[fl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g3  =  Plotffl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g4  =  Plot[fl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

g5  =  Plot[fl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g6  =  Plotffl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g7  =  Plotfff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g8  =  Plotfff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g9  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=  1  \[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

glO  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=l\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 
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gl  1  —  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Dram  (R=50cm,rpm=8,D=  1  \[Mu])" , 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  .0000001; 

DD  =  2; 
gl2  =  Plot[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  (thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Dram  (R=50cm,rpm=0,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g22  =  Plot[fl,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Dram  (R=50cm,rpm=l,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g32  =  Plot[fl,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g42  =  Plot[fl,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

g52  =  Plotffl,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Dram  (R=50cm,rpm=6,D==2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g62  =  Plot[fl,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g72  =  Plotfff,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g82  =  Plotfff,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g92  =  Plot[ff,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl02  =  Plot[ff,  (thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=2\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 
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gl  12  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=2\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  .0000001; 

DD  =  3; 
gl3  =  Plot[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  {thour,  0.001,  50}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g23  =  Plotffl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g33  =  Plotffl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g43  =  Plot[fl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

g53  =  Piotffl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g63  =  Plot[fl,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g73  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g83  =  Plotfff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g93  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=3\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl03  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=3\[Mu])M, 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 
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gl  1 3  =  Plot[ff,  {thour,  0.001,  50}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=3\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  .0000001; 

DD  =  4; 
gl4  =  Plot[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  f2]],  (thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm:=0,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g24  =  Plotffl,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g34  =  Plotffl,  (thour,  0.001, 25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm==2,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g44  -  Plotffl,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6)‘ 

g54  =  Plotffl,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g64  =  Plotffl,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm::=8,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g74  =  Plotfff,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g84  =  Plotfff,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g94  =  Plotfff,  (thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl04  =  Plotfff,  (thour,  0.001,  25}, 

.  PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=4\[Mu])", 
AxesLabel  ->  ("thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 
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gl  14  =  Plot[ff,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=4\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  .0000001; 

DD  =  5; 
gl  5  =  Plot[ 

If[h/r  <=  hrlow,  fl,  If[hrhigh  <=  h/r,  f3,  £2]],  {thour,  0.001,  25}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g25  =  Plot[fl,  {thour,  0.001, 25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g35  =  Plot[fl,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g45  =  Plot[fl,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R:=50cm,rpm=4,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

g55  =  Plotffl,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g65  =  Plotffl,  {thour,  0.001, 25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g75  =  Plot[ff,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g85  =  Plot[ff,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g95  =  Plot[ff,  {thour,  0.001, 25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl05  =  Plotfff,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 
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gl  15  =  Plot[ff,  {thour,  0.001,  25}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=5\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  .0000001; 

DD  =  6; 

gl6  =  Plot[If[h/r  <=  hrlow,  fl,  iprhigh  <=  h/r,  O,  f2]],  {thour,  0.001 ,  7}, 
PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=0,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  1; 

g26  =  Plot[fl,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=  1  ,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g36  =  Plot[fl,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D:=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g46  =  Plot[fl,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

g56  =  Plotffl,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

g66  =  Plot[fl,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=6\[Mu])M, 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm=  1; 

g76  =  Plot[ff,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=l,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  2; 

g86  =  Plot[ff,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=2,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  4; 

g96  =  Plotfff,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=4,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  6; 

gl06  =  Plot[ff,  {thour,  0.001,  7}, 

PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=6,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
rpm  =  8; 

gl  16  =  Plotfff,  {thour,  0.001,  7}, 
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PlotLabel  ->  "Rotating  Drum  (R=50cm,rpm=8,D=6\[Mu])", 
AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Plot[KK,  {DD,  1,  10}, 

PlotLabel  ->  "  Electrostatic  Dispersion(Ncharge=10*DA2)", 
AxesLabel  ->  ("D(\[Mu])",  "Ke  PlotRange  ->  {0,  10A-6}, 
PlotPoints  ->  20]; 

Show[gl,  g2,  g3,  g4,  g5,  g6,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Moss  Theory  (rpm=0,l,2,4,6,8;D=l\[Mu])"]; 
Show[g7,  g8,  g9,  glO,  gll,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Ncharge=T0*DA2  (rpm=l,2,4,6,8;D=l\[Mu])"]; 
Show[gl2,  g22,  g32,  g42,  g52,  g62,  PlotRange  ->  (0,  1}, 

PlotLabel  ->  "  Moss  Theory  (rpm=0,l,2,4,6,8;D=2\[Mu])"]; 
Show[g72,  g82,  g92,  gl02,  gl  12,  PlotRange  ->  {0,  1), 

PlotLabel  ->  "  Ncharge=10*DA2  (rpm=l,2,4,6,8;D=2\[Mu])"]; 
Show[gl3,  g23,  g33,  g43,  g53,  g63,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Moss  Theory  (rpm=0,l,2,4,6,8;D=3\[Mu])"]; 
Show[g73,  g83,  g93,  gl03,  gll3,  PlotRange  ->  (0,  1}, 

PlotLabel  ->  "  Ncharge=10*DA2  (rpm=l,2,4,6,8;D=3\[Mu])"]; 
Show[gl4,  g24,  g34,  g44,  g54,  g64,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Moss  Theory  (rpnr=0, 1 ,2,4,6,8;D=4\[Mu])"] ; 
Show[g74,  g84,  g94,  gl04,  gl  14,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Ncharge=10*DA2  (rpm=l,2,4,6,8;D=4\[Mu])"]; 
Show[gl5,  g25,  g35,  g45,  g55,  g65,  PlotRange  ->  (0,  1}, 

PlotLabel  ->  "  Moss  Theory  (rpm=0,l,2,4,6,8;D=5\[Mu])"]; 
Show[g75,  g85,  g95,  gl05,  gl  15,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Ncharge=10*DA2  (rpm=l,2,4,6,8;D=5\[Mu])"]; 
Show[gl6,  g26,  g36,  g46,  g56,  g66,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Moss  Theory  (rpm=0,l,2,4,6,8;D=6\[Mu])"]; 
Show[g76,  g86,  g96,  gl06,  gl  16,  PlotRange  ->  {0,  1}, 

PlotLabel  ->  "  Ncharge=10*DA2  (rpm-l,2,4,6,8;D=6\[Mu])"]; 

MATHEMATICA  file  RadiusOfGyration.nb 

(*using  equivalent  aerodynamic  diameter  so  density  =  1  *) 
rho  =  1 ; 

mu=  1.83*1 0A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*. 07 1  *  1 0A-4/d)( 1.257  +  ,4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
p  =  Exp[-2*tau*t*omegaA2]; 
rO  =  tau  980/omega; 
r  =  50; 

gl  =  Plot3D[r0,  (DD,  .001,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Radius  of  Gyration  r0(cm)", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "r0"},  PlotPoints  ->  30]; 
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Show[%,  PlotRange  ->  {0,  3}]; 

volfractleft  =  If[rO  >  r,  0,  (r  -  rO)A2/rA2]; 

gl  =  Plot3D [volfractleft,  {DD,  .001,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "N/NO  After  One  Revolution;  R=50cm", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/N0  "},  PlotPoints  ->  30]; 
Show[%,  PlotRange  ->  {0.85,  1}]; 

MATHEMATICA  file  RatioKNEOoverKeqO.nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*.071*10A-4/d)(1.257  +  .4*Exp[-l.l*d/(2*.071*10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
r  =  50; 

nn  =  l/(Exp[2*tau*t*omegaA2]* 

(1/nnO  +  KK/(2*tau*omegaA2))  -  KK/(2*tau*omegaA2)); 

(*KK  =  8*Pi*tau*qpart*qcloudavg/3/mm; 
mm  =  rho*Pi*dA3/6; 
qpart  =  ncharge*4.8*10A-10; 
ncharge  =  10*DDA2; 
qcloudavg  =  10*DDsquaredmean; 

DDsquaredmean  =  9;*) 
ff  =  nn/nnO; 

ffoverffO  =  ff/Exp[-2  *  tau*  t*omegaA2] ; 
thour  =  48; 

KK  =  0; 
nnO  =  1000; 

gl  =  Plot3D[ff,  {DD,  1,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Kcoag  +  Kelec  =  0  ;(No=1000,t=48hr)", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 

KK  =  10A-9; 

gl  =  Plot3D[ff,  {DD,  1,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Kcoag  +  Kelec  =  10A-9  ;(No=1000,t=48hr)", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"},  PlotPoints  ->  20]; 
gl  =  Plot3D[ffoverff0,  {DD,  1,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Nl(Kcoag+Kelec=10A-9)/N2(K=0);  No=1000,t=48hr", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N1/N2  "},  PlotPoints  ->  20]; 
nnO  =  500; 

gl  =  Plot3 D [ffoverffO,  {DD,  1,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Nl(Kcoag+Kelec=10A-9)/N2(K=0);  No=500,t=48hr", 
AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "N1/N2  "},  PlotPoints  ->  20]; 
nnO  =  1000; 

KK  =  10A-8; 

gl  =  Plot3D[ff,  {DD,  1,10},  {rpm,  1,  8}, 
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PlotLabel  ->  "Kcoag  +  Kelec  =  10A-8  ;(No=1000,t=48hr)", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N/No"} ,  PlotPoints  ->  20]; 
gl  =  Plot3 D [ffoverffO,  {DD,  1,  10},  (rpm,  1,  8}, 

PlotLabel  ->  "N 1  (Kcoag+Kelec=  1 0A-8)/N2(K=0);  No=1000,t=48hr", 
AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "N1/N2  "},  PlotPoints  ->  20]; 
nnO  =  500; 

gl  =  Plot3D[ffoverff0,  {DD,  1,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Nl(Kcoag+Kelec=10A-8)/N2(K=0);  No=500,t=48hr", 
AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "N1/N2  "},  PlotPoints  ->  20]; 
nnO  =  700; 

KK=  10A-9; 

gl  =  Plot3D[ffoverff0,  {DD,  1,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Nl/N2\[Congruent]N(K=10A-9)/N(K=0);  No=700,t=48hr", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N1/N2  "},  PlotRange  ->  {.85,  1}, 
PlotPoints  ->  20]; 
nnO  =  700; 

KK=  10A-8; 

gl  =  Plot3D[ffoverff0,  {DD,  1,  10},  {rpm,  1,8}, 

PlotLabel  ->  "Nl/N2\[Congruent]N(K=10A-8)/N(K=0);  No=700,t=48hr", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N1/N2  "},  PlotRange  ->  {.4,  1}, 
PlotPoints  ->  20]; 
thour  =  24; 
nnO  =  700; 

KK  =  8*Pi*tau*qpart*qcloudavg/3/mm; 
mm  =  rho*Pi*dA3/6; 
qpart  =  ncharge*4.8*10A-10; 
ncharge  =  10*DDA2; 

qcloudavg  =  10*DDsquaredmean*4.8*10A-10;  (^depends  on  dist  funct*) 
DDsquaredmean  =  9; 

gl  -  Plot3D[ffoverff0,  {DD,  1,  10},  {rpm,  1,  8},  PlotLabel  ->  "\!\(\* 
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StyleBox[\"  A"  ,\n\"  SmallT  ext\"  ,\nF  ontSize->  1 0]\)\!\(\* 
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AxesLabel  ->  {"D(\[Mu])",  "rpm",  "N1/N2  "},  PlotPoints  ->  20]; 
gl  =  Plot3D[KK,  {DD,  \,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Kelec  withNcharge:=10*DA2;  No=700,t=24hr", 

AxesLabel  ->  {"D(\[Mu])",  "rpm",  "K  "},  PlotPoints  ->  20]; 

MATHEMATICA  file  DiffusionDeposition.nb 

t  =  thour*3600; 
r  =  50; 

f  =  (6/PiA2)Sum[Exp[-DiffcoePPiA2*iA2*t/rA2]/iA2,  {i,  1,99}]; 

Diffcoef  =  1 0A- 1 ; 

gl  =  Plot[f,  {thour,  0.001,  50}, 

PlotLabel  ->  "Diffusion  Solution  (R=50cm,DC=10A-lcmA2/s)", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Diffcoef  =  10A-2; 

g2  =  Plot[f,  {thour,  0.001,  50}, 

PlotLabel  ->  "Diffusion  Solution  (R=50cm,DC=10A-2cmA2/s)", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Diffcoef  =  10A-3; 

g3  =  Plot[f,  {thour,  0.001,  50}, 

PlotLabel  ->  "Diffusion  Solution  (R=50cm,DC=T0A-3cmA2/s)", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Diffcoef  =  10A-4; 

g4  =  Plot[f,  {thour,  0.001,  50}, 

PlotLabel  ->  "Diffusion  Solution  (R=50cm,DC=10A-4cmA2/s)", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 

Diffcoef  =  10A-5; 

g5  =  Plot[f,  {thour,  0.001,  50}, 

PlotLabel  ->  "Diffusion  Solution  (R=50cm,DC=10A-5cmA2/s)", 

AxesLabel  ->  {"thour",  "N/No"},  PlotPoints  ->  20]; 
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Show[gl,  g2,  g3,  g4,  g5,  PlotRange  ->  {{0,  50},  (0,  1}},  PlotLabel  ->  "\!\(\* 
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MATHEMATICA  file  DiffusionDepositionSurf.nb 

rho  =  1 ; 

mu  =  1.83*10A-4; 
d  =  DD*10A-4; 
omega  =  2*Pi*rpm/60; 
t  =  thour*3600; 

CC  =  1  +  (2*  .07 1  *  1 0A-4/d)( 1.257  +  .4*Exp[- 1 . 1  *d/(2*.07 1  *  10A-4)]); 
tau  =  rho*dA2  CC/(18*mu); 
r  =  50; 
b  =  r; 

vcent  =  tau*omegaA2*r; 
f  =  (1  -  Exp[-vcent*b/Diffcoefj); 

DD  =  2; 
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rpm  =  4; 

Diffcoef  =  1; 

gl  -  Plot3D[f,  {DD,  0.001,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Deposition  without/with  diffusion;  D=lcmA2/s", 
AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "WO/W"},  PlotPoints  ->  20]; 

Diffcoef =  10A-1; 

gl  =  Plot3D[f,  {DD,  0.001,  10},  {rpm,  1,  8}, 

PlotLabel  ->  "Deposition  without/with  diffusion;  D=.lcmA2/s", 
AxesLabel  ->  {"D(\[Mu])",  "rpm",  "WO/W"},  PlotPoints  ->  20]; 

Diffcoef  =  10A-2; 

gl  =  Plot3D[f,  {DD,  0.001,  10},  {rpm,  1,8}, 

PlotLabel  ->  "Deposition  without/with  diffusion;  D=.01cmA2/s", 
AxesLabel  ->  {"D(\[Mu])M,  "rpm",  "WO/W"},  PlotPoints  ->  20]; 
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